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Recent studies of the Q 2 dependence of ep scattering in the large x region and in the Q 2 
range: 1 < Q 2 < 30 Gev 2 , confirm the validity of the phenomenon of quark-hadron duality - 
the similarity of the deep inelastic (parton) and resonance (hadron) spectra - for values of the 
invariant mass, W 2 > 2.4 GeV 2 . At lower values of W 2 , duality is found to be significantly 
violated by an amount that cannot be parametrized in terms of the first few terms of a series 
of power corrections. We present a dynamical model that explains the Q 2 dependence of 
the data: at low W 2 , non-partonic components given by color neutral clusters dominate the 
cross section and the Q 2 dependence is governed by their mass spectrum, predicted within 
the preconfiment property of QCD; at large W 2 the structure function is determined by a 
convolution of the cluster mass spectrum with the parton distributions. 



1 Introduction 



Deep inelastic scattering (DIS) experiments provide a wealth of information on the structure of 
hadrons, from a detailed understanding of their partonic structure, to precision measurements of 
as- In order to reveal the partonic structure of hadrons, the four- momentum transfer squared, 
Q 2 , must be large, namely well above 1 GeV 2 . As Q 2 is lowered partonic components become less 
likely to be observed, and a transition to a regime eventually dominated by non-partonic degrees 
of freedom, and therefore outside the range of applicability of pQCD, is expected to occur. The 
onset of this regime has been observed in a large number of experiments on ep scattering in 
the low Bjorken x and large invariant mass, W 2 , region, where it is a well accepted fact that 
non-partonic components should replace the partonic structure at low Q 2 . The data indeed 
support a sharp transition between the partonic and non-partonic regimes at Q 2 ~ 1 GeV 2 , 
for 10 -5 < x < 10~ 3 . At large x the proton structure function is clearly dominated by non- 
partonic components - the nucleon resonances - up to relatively large Q 2 (Q 2 < 20 GeV 2 ). An 



equivalence was singled out, however, between the average of the resonance spectrum, written in 
terms of Mellin moments, and the DIS structure function measured at larger Q 2 . The moments 
in the two regimes were found to differ by perturbative corrections and relatively small power 
corrections!. It was conjectured that a form of duality resulted from the cancellation of higher 
order terms in the twist expansion that would otherwise be expected to dominate the cross 
section at x — > 1, or as more exclusive states are produced. This view has been since considered 
the most natural one. From a slightly different perspective, a series of recent papers has been 
devoted to local auark-hadron duality and its violations in semi-leptonic decays, and r decays 
(for a review seea). In particular, it was shown incl that local duality violations can be traced 
to the asymptotic nature of the operator product expansion (OPE), namely, to the behavior of 
operators of both higher dimension and higher twist. This has lead the way to specific models 
based on instantons, and on large N c QCD in (1 + 1) dimensions. Approaches consistent withEra 
could shed some light on duality in DIS, where quark-hadron duality has been explored so far 
within the context of quark models 0, thus avoiding the basic questions related to the dynamical 
nature of QCD. In this contribution, we discuss a possible avenue: Our starting point is similar 
to Bel, in that the background of our model is the OPE within which we pursue a connection 
between the resonance region and the higher twist operators. Crucial for the construction of 
our model is an accurate analysis of recent dataQ conducted inQ. In the following sections we 
summarize the results found inQ, we introduce our model, and we present some initial results. 



2 Perturbative QCD Analysis of Parton-Hadron Duality 

In@ the Q 2 dependence in the resonance region was extracted by first considering the average 
of the resonances over £ = 2x/(l + \/l + 4M 2 x' 2 / Q' 2 ) , which properly takes into account Target 
Mass Corrections (TMC). The averaging procedure, described in detail in i, yields a smooth 
curve in £ which fits the resonance data with a x 2 /d.o.f. between 0.8 and 1.1. It represents an 
alternative analysis to the ones using moments. The fit was performed in bins of W 2 , centered 
at: Wft = 1.6, 2.3, 2.8, 3.4 GeV 2 , respectively. Our subsequent study was aimed at establishing 
whether it is possible to define a breakpoint where jpQCD no longer applies and a transition 
occurs, similar to what observed in the low x regimeQ. The analysis involved a number of steps 
similar to recent extractions of power corrections from inclusive dataBot3, namely the form 

F?»(x, Q 2 ) = F^ CD+TMC (x, Q 2 ) + H[X n f ] + 0(1/Q 4 ), (1) 

was adopted, where F2® CD+ ™ C (x, Q 2 ) is the twist- 2 contribution, including the kinematical 
TMC; the other terms in the formula are the dynamical power corrections, formally arising 
from higher order terms in the twist expansion. Both jtpQCD+tmc and H 

were extracted from 

the data at large x. The shape of the initial NS PDFs was found to be well constrained at 
variance with the singlet and gluon distributions at low Q 2 , whose shape is strongly correlated 
with the value of as- NNLO corrections were not included. These would introduce further 
theoretical uncertainties. In fact, the question of whether they can "mimick" the contributions 
of higher twists, including the uncertainties due to the well known scale/scheme dependence of 
calculations, within the current precision of data is still a subject of intense investigations!!. Large 
x resummation was performed directly in x space by replacing the Q 2 scale with a z-dependent 
one, W 2 = Q 2 (l— z)/z0. It was found that over the range 0.45 < x < 0.85, large x resummation, 
and TMC improve the agreement with the data. We parametrized the remaining discrepancies 
through the Higher Twist (HT) coefficient H{x,Q 2 ) = F$ QCD+TMC (x, Q 2 )C H t{x). In Fig. la 
we show Cht-i extracted from: DIS data with W 2 > 4 GeV 2 , from the resonance region (W 2 < 4 
GeV 2 ), and over the entire range of W 2 . The figure also shows the range of extractions previous 
to the current oneQ. While the large W 2 data track a curve that is consistent with the 1/W 2 



behavior expected from most models, the low W 2 data yield a much smaller value for Cht an d 
they show a bend-over of the slope vs. x. This surprising effect is not a consequence of the 
interplay of higher order corrections and the HT terms, but just of the extension of our detailed 
pQCD analysis to the large x, low W 2 kinematical region. In other words, we unraveled a 
Q 2 dependence that seems to deviate from the common wisdom developed since the pioneering 
analysis ofi, or, in the language of @, we observe a violation of global duality. Similar results, 
expressed through a Q 2 dependent HT coefficient, were also obtained intil. 



3 Large N c Model for Initial Parton Evolution 

We propose a simple dynamical model for the structure function in the low W 2 (W 2 < 4 GeV 2 ) 
and low Q 2 (Q 2 < 10 GeV 2 ) regime, where non-partonic configurations are expected to be 
dominant. The DIS cross section for ep scattering is proportional to the structure functions 
F^ 2 ^(x,Q 2 ) which in turn, measure combinations of the parton longitudinal momentum dis- 
tributions, qi(x,Q 2 ), i = u,d,s,... at the scale Q 2 . In the standard approach to DIS the Q 2 
dependence of F2 is described by the pQCD evolution equations El, whose numerical solution 
requires parametrizing the input distributions at an initial scale Q 2 where pQCD is believed 
to be still applicable. Q 2 serves as a boundary between the perturbative and non-perturbative 
domains, although its value is somewhat arbitrary (in current parametrizations Q 2 ~ 0.4 — 10 
GeV 2 ). We refer to this situation as the "fixed initial scale" description, and we write explicitely 
the dependence of qi(x, Q 2 , Q 2 ) on Q 2 . A simple kinematical argument shows that Q 2 is related 
to the invariant mass squared of the proton remnant after a parton is emitted, by: M\ ~ Q 2 /x. 
At large x M\ ~ Q 2 (at low x this simple kinematical observation, and the factorization prop- 
erties of the diffractive part of F2, support the idea that the parton is emitted from a large mass 
object, indentified with a soft pomeron 0). We, therefore, explore the possibility that partons 
are not emitted directly from the nucleon, but that, before the pQCD radiative processes are 
initiated, a semi-hard phase occurs where the dominant degrees of freedom are color neutral 
clusters with a mass distribution peaked at: /u 2 eafc « Q 2 . As a result, the nucleon structure 
function is related to the quark distribution by a smearing of the initial Q 2 , namely 

r W 2 fj.,2 

F 2 (x,Q 2 ) = xJ2e 2 / ^P^ 2 ) qi (x,Q 2 ^ 2 ), (2) 

where P(n 2 ) {Ppeaki^ 2 ) ~ P(Qo)) 1S the clusters' mass distribution, and the sum is extended 
to valence quarks only since we are describing the large x region. Eq.(^) expresses the fact that 
the initial stage of pQCD evolution is characterized by color neutral clusters of variable mass, 
from which the hard scattering parton will emerge, in a subsequent stage of the interaction. 

Our model is derived within the framework of the large N c approximation, in analogy to the 
cluster hadronization schemes implemented in the HERWIG Monte Carlo simulation Ej. In this 
scheme, hadronization proceeds as prescribed by the pQCD property of preconfinement of color 
@: at the end of the parton's pQCD evolution, color singlets axe formed with a Q 2 - independent 
mass (and spatial) distribution. In practical implementations^^, all gluons left at the hadroniza- 
tion scale, are "forcibly" , or non-perturbatively, split into qq pairs. It is this modification of the 
evolution equations that allows for the local parton-hadron conversion through preconfinement 
of colour: each color line "color-connects" e.g. a quark to an anti-quark, forming a color singlet. 
The color-singlet clusters are then fragmented into hadrons. In DIS the transiton hadrons — > 
quarks — > hadrons, is complicated both by initial state radiation and by the presence of the beam 
cluster formed from the remnant of the initial hadron. This produces an additional rescattering 
term in the function P in Eq. (2) 113. 

As a preliminary study, we considered both the low and the very large W 2 limits of Eq.@. At 
low W 2 , F2(W 2 , Q 2 ) « P(Q 2 ), namely it is described by the behavior of the cluster distribution 
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Figure 1: (a) Higher Twist coefficient from Eq.(2); (b) Q 2 dependence at fixed W 2 = 1.6 GeV 2 . 

function. At very large W 2 , P pea k determines the value of the initial Q 2 . In Fig. 1(b) we present 
our result for F2 at W 2 = 1.6 GeV 2 , where the cluster distribution was obtained from HERWIG 
113 (triangles), whereas the red curve is an analytic calculation of the Sudakov-type behavior 
valid at large Q 2 . 
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